Chiral Disorder and Random Matrix Theory with Magnetism 
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We revisit the concept of chiral disorder in QCD in the presence of a QED magnetic field \eH\. 
Weak magnetism corresponds to \eH\ < 1/p 2 with p ~ 1/3 fm the vacuum instanton size, while 
strong magnetism the reverse. Asymptotics (ultra-strong magnetism) is in the realm of perturbative 
QCD. We analyze weak magnetism using the concept of the quark return probability in the diffusive 
regime of chiral disorder. The result is in agreement with expectations from chiral perturbation 
theory. We analyze strong and ultra-strong magnetism in the ergodic regime using random matrix 
theory including the effects of finite temperature. The strong magnetism results are in agreement 
with the currently reported lattice data in the presence of a small shift of the Polyakov line. The 
ultra-strong magnetism results are consistent with expectations from perturbative QCD. We suggest 
a chiral random matrix effective action with matter and magnetism to analyze the QCD phase 
diagram near the critical points under the influence of magnetism. 



INTRODUCTION 

Chiral quarks in the QCD vacuum break spontaneously 
chiral symmetry. Key contributors to this spontaneous 
breaking are instanton and antiinstanton fluctuations 
with left handed and right handed zero modes attached 
to them [IHS] (and references therein). The random na- 
ture of the instanton and antiinstanton fluctuations cause 
these zero modes to spread near zero virtuality. This fun- 
damental property of the QCD quark spectrum is cap- 
tured by the Banks-Casher relation |4 



7r/9 4 (0) 



(1) 



that ties the quark condensate to the quark spectral den- 
sity p 4 (A) at zero virtuality. Note the positive convention 
for the chiral condensate. In the QCD vacuum with in- 
stanton density 114 and size p, the quark condensate is 
expected to scale as £4 « 114/5. Using the Gell-Mann 
Oakes Renner relation m 2 F 2 = 2to£4, (fil) is just the 
analogue of the Kubo formula for the DC conductivity 
in metals 



cr c = F'11-k = D Pi (0) 



(2) 



The pion decay constant F n defines the chiral conduc- 
tivity, with the chiral diffusion constant D = F^/Y.4. 
Through ( l]J2 ) chiral quarks trapped in an Euclidean 4- 
volume V4 = L 4 behave much like electrons in disordered 
metallic grains with At = D/L 2 playing the role of the 
Thouless energy [SHE]- Quark with virtualities X < Xt 
are in the ergodic regime and are well described by ran- 
dom matrix theory, while quarks with virtualities X > Xt 
are in the diffusive regime which is amenable to chiral 
perturbation theory. In the presence of a fixed external 



magnetic field H the chiral disorder is affected. In this 
paper we address the details of these changes. 

Weak magnetism is mainly affecting the diffusive prop- 
erties of chiral fermions through a renormalization of 
the low energy parameters, e.g. £4,-^,771^. Although 
the chiral expansion in the massless limit relies on 
\eH\/(ATrF 1T ) 2 < 1 and therefore suggests that weak mag- 
netism operates in the realm of eH\ < 1 GeV 2 , the valid- 
ity range is substantially smaller. Indeed, since instan- 
tons may be at the origin of the spontaneous breaking 
of chiral symmetry we expect weak magnetism to break 
down when the instanton size p is resolved. As a result, 
the range of weak magnetism is \eH\ < 1/p 2 « 1/3 GeV 2 
and will be referred to as diffusive. 

Strong magnetism corresponds to \eH\ > 1/p 2 ~ 
1/3 GeV 2 with a magnetic field starting to dwarf the 
chromomagnetic field carried by the instanton. This 
range will be referred to as ergodic. As a result, light 
and chiral quarks in the QCD vacuum switch from lock- 
ing their spin to the instanton color to lining their spin 
with the external magnetic field. They form Landau or- 
bits which are gapped by \eH\ near zero virtuality. The 
exception is the lowest Landau level (or LLL for short), 
where the magnetic contribution cancels the zero point 
motion. As a result, the LLL with its huge degeneracy 
piles up at zero virtuality becoming a potential contribu- 
tor to the spontaneous breaking of chiral symmetry. The 
LLL is inherently 2-dimensional. This phenomenon of di- 
mensional reduction from 4 to 2 in Euclidean space with 
accumulation near zero virtuality of the LLL favors the 
spontaneous breaking of chiral symmetry through any 
residual interaction. This phenomenon is known as mag- 
netic catalysis [5] (and references therein). For instance, 
residual longitudinal interactions in the LLL through the 
instantons will cause the spontaneous breaking of chiral 
symmetry, 



S 4 = |e£T|E 2 » \eH\yftup (3) 

with typically n 4 ss 1/fm the vacuum instanton density. 
Ultra-strong magnetism is solely characterized by Lan- 
dau levels with huge degeneracy and small sizes. As a 
result, many Landau orbits can fit in a single instanton 
making the concept of instanton zero modes not particu- 
larly useful. In many ways, this picture resembles that of 
the initial color glass condensate with many wee partons 
saturating the transverse size of a colliding nucleus [TU] 
(and references therein). Perturbation theory becomes 
the rule and will cause the LLL to spread near zero vir- 
tually. The spontaneous breaking of chiral symmetry 
through perturbative gluons follows with 



E 4 = | ej ff|E 2 « |eiJ| 3 / 2 



(4) 



as expected from dimensional arguments. We expect this 
scaling to take place for \eH\ > 10 / p 2 « 3 GeV 2 or when 
about 10 Landau orbits can fit within the instanton trans- 
verse size thereby dwarfing the zero modes. 

In section 2, we discuss weak magnetism in the diffusive 
regime of 4-dimensional chiral quarks trapped in U 4 using 
the concept of the quark return probability. The result is 
in line with the one from chiral perturbation theory. In 
sections 3 and 4, we discuss ultra-strong and strong mag- 
netism in the ergodic regime where the instanton size is 
resolved through the use of random matrix theory. The 
results arc in line with the expectations of dimensional 
reduction from the LLL. In section 5 we compare our re- 
sults to current lattice simulations |llj . In particular we 
show that a small shift in the trivial Polyakov holonomy 
is magnified by the chiral transition and may account 
for the anti-catalysis observed on the lattice as recently 
suggested in |12) . In section 6, we show how random ma- 
trix models for QCD with magnetism can be constructed 
to shed light on the comparison with widely used con- 
stituent quark models of the NJL type, and help analyze 
a number of issues in the QCD phase diagram with mag- 
netism. Our conclusions and prospects are in section 7. 



DIFFUSIVE REGIME WITH WEAK 
MAGNETISM 

The dynamics of QCD light quarks in a 4-dimensional 
Euclidean box is chiral and diffusive in the long wave- 
length limit. The diffusive nature is best captured by 
the quark return probability [5] 



P(0,r) = (\u + (t)u(0) +d+(r)d(0)| 



(5) 



for 2 light flavors. In the absence of magnetism, the vac- 
uum is isospin symmetric and ([5]) is dominated by the 
the triplet of pions in an Euclidean box 



P(0,r)=2(P (0,r)+P±(0,r)) 



with 



Po,±(0,r) 



yV-DQVi 



(6) 



(7) 



for a triplet of massless pions. The sum is over the 
isotriplet of charged pions or diffuson modes with mo- 



menta Q^ 



^tt/L in a periodic V4 — L 4 Euclidean 



box. The quark return probability ([5]) in the chiral limit 
obeys a sum rule 



1 r°° 

E 4 = - lim lim — / P(0,r)dr 

rn— >Q Vi— >oo V4 Jq 



(8) 



In the presence of magnetism, the vacuum is no longer 
isospin symmetric. As a result the free chargeless pion 
7T° remains gapless, while the free charged pions ir^ are 
gapped in Landau levels (LL). For a constant magnetic 
field with A 1 = —Hx2, the LL are 



X 2 ± = \eH\(2n + l) 



2 1 2 
P3+P4, 



(9) 



with degeneracy <fr — eHL 2 /2ir. The change in the quark 
return probability follows from the change in the diffuson 
propagator for the charged pion modes 



^ eHL 2 



n,m,k 



2tt 



-_Dr[(n+l/2)2eff+(m 2 +fc 2 )(27r/L) 2 +m|] 



(10) 



The change in the quark return probability is the change 
in the charged diffuson modes and is captured by the 
difference 



/--oc 

/= / [P±(ff,r)-P±(0,T)]dT 

Jo 



(11) 



In the chiral limit, replacing the sums over free momenta 
by integrals and summing explicitly over the Landau lev- 
els of the charged pions (i.e. the diffusons) we obtain 



1 



eHV 



16n 2 D J \zsmh 



72 ,,/: 



eHV 
16n 2 D 



In 2 



(12) 



Using the value of the diffusion constant we arrive at 



E 4 (ff) - £ 4 (0) 1 



eH In 2 
16tt 2 P2 



(13) 



which is the result of chiral perturbation theory in leading 
order [13] . As expected, the diffusive regime is the regime 
of chiral perturbation theory albeit in a finite Euclidean 
box [14] . Chiral perturbation theory in a finite box was 
formulated systematically in [13] , 



ERGODIC REGIME WITH ULTRA-STRONG 
MAGNETISM 

For strong magnetism the chiral disorder in 4 dimen- 
sions with an isotriplet of pions transmute to a chiral 
disorder in quasi- 2-dimensions with a chargless pion, as 
the charged diffusons become gapped. The transmuta- 
tion takes place when the magnetic field resolves the in- 
stanton size or \eH\ > l/p 2 ~ 1/3 GeV 2 . In this regime, 
magnetism dwarfs the instanton chromo-magnetism with 
Landau levels becoming the lore. The LLL dominates the 
infrared physics near zero virtuality triggering the catal- 
ysis of chiral symmetry breaking by magnetism. 

Indeed, for strong magnetism the free quark spectrum 
in 4-dimcnsions is given by Landau levels 



\± a (w, k z ) = ±^/uj 2 + k 2 z + \eH\(2n + 1 - s) + m 2 

(14) 
with current quark masses. For simplicity we will con- 
sider Np — 1 with charge e unless specified otherwise. 
Each (u,k z ) branch is 



^H\L x L y /2ir=\H\L x L y /(hc/e) 



(15) 



degenerate. ( 15 ) measures the magnetic flux in units 



of the quantum flux he/ e and is integer valued on a 
4-dimensional lattice with periodic (toroidal) boundary 
conditions. The Lowest Landau branches n = 0, s = 1 



\± 1 (LJ,k z ) = ±y/bJ 2 + k* + 



rn- 



(16) 



are paired by chirality and dominate in the infrared limit. 
They are 2(j> degenerate. The higher Landau levels are 
gapped at zero virtuality. We note that at finite temper- 



ature the w-spectrum in ( 16 1 is dominated by the lowest 
Matsubara mode u — > ujq = ttT (16j . This approximation 
will be used to probe critical points in phase diagrams 
with magnetism. 

While the color interactions are involved within and 
between Landau levels, a model independent understand- 
ing can be achieved through the Banks-Casher formulae 
used earlier. In the presence of strong magnetism the 
localization of the quark spectrum in the transverse hy- 
perplane to the magnetic field suggests 



(*ltf 4 ) = 7rp 4 (0) = y 



2<\> 
4 AA 



\eH\ 



f3L z A\ 



(17) 



where AA is the typical level spacing near zero virtuality. 
In a spontaneously broken QCD vacuum with no mag- 
netism AA « I/V4. The last relation in (17) reflects on 



the dimensional reduction of the quark spectral function 
near zero virtuality from 4 to 2 dimensions, 



P4 (0) = \eH\p 2 (0) 
and therefore the condensate relation 



(18) 



^i* 4 ) = \eH\ (^2) (19) 

In the infrared regime, this suggests the field redefinition 



* 4 (a; , x, y, z)) -> s/\eH\ # 2 (z , z) 



(20) 



with the pertinent canonical dimensions made explicit. 

The breaking of chiral symmetry in terms of ^ 2 is 
quasi-2-dimensional. It does not upset the Mermin- 
Wagner theorem. Indeed, the random gluonic inter- 
actions that cause the spontaneous breaking of chiral 
symmetry and therefore the redistribution of the 2- 
dimensional quark spectrum near zero virtuality are still 
blind to magnetism and fully 4-dimensional. 

In quasi- 2-dimensions the free quark fields ^2 are dom- 



inated by the branches ( 16 ) and any residual gluonic 
interaction will cause them to redistribute. For ultra- 
strong magnetism with 

\eH\> 10/p 2 «3GeV 2 

the interactions among the quarks in the LLL are mainly 
perturbative as the instanton quasi-zero modes become 
obsolete. Assuming the residual interactions in the in- 
frared for the lowest landau level (LLL) and lowest Mat- 
subara frequency to be totally random but chiral, the 
partition function for these modes is best captured by a 
new chiral RMM 



Zlll(™,T, 4>) = 



771 IUJ{) + iA 



(21) 



where A is a <f> x <fi real matrix with a Gaussian distri- 
bution. Recall that we have set Np — 1. The novelty 



of (21) is in the scaling of the Gaussian distribution of 
the matrix elements with magnetism. It captures the 
random interactions induced by plain perturbative glu- 
ons between the (f> degenerate quarks S&2 in the LLL and 
amounts to the 4-fermi interaction 



\eH\ 



*|*2 



as expected from dimensional reduction. The Gaussian 
weight enforces a level spacing between the eigenvalues 
of the chiral matrix to be of order 



AA 



y/W\ 



(22) 



By the Banks-Casher formulae we expect the chiral con- 
densate in a symmetric box V4 = L 4 



IT 2(f) 



1 



*W=™<°> = ££~^'" ! 



(23) 



in agreement with the dimensional arguments presented 
in Q. This result is born out by explicit calculations as 
we now show. 



The quark condensate follows through ( 19 ) as 



*I* 4 > = lim 



\eH\ dWL LLL 



i->o /3L Z dm 



(24) 



In the large (j> limit, the partition function in (24) is an- 



alyzed through standard arguments by fermionizing the 



determinant in (21), performing the Gaussian integral 



and then bosonizing the resulting 4-fermi interaction us- 
ing the Hubbard-Stratonovitch transform, 



^LLL 



(m, T, <j)) = 



dbe-T^ hth det« 



b + m iujQ 
iujQ b^ + m 



(25) 



Inserting (25) in (24) yields 



*4*4) = lim lim 



\eH\ 



m— >0L Z — >oo pL z 



tr 



b + m iujQ 
iujQ b^ + m 



(26) 



The averaging is carried using the measure in (25). In 
the saddle point approximation, the result is 



*I*4 



\eH\ 20 



which can be re-written as 



l-u%/\eH\ (27) 



are consistent with dimensional arguments for asymp- 
totically strong magnetism, and the perturbative QCD 
arguments presented in [131 . This analysis ignores the 
back-reaction of the quarks on the gauge configurations 
in the saddle point approximation. The back-reaction 
maybe important as we detail below. 

Finally, we note that a more general RMM for all Lan- 
dau levels can be written similarly by insisting on the 
free Landau level spectrum for A = 0. Specifically 



/oo 
dAe -W/|eH|)1v(AtA)-Q JJ 
CI r, L 1 



xdet 



n=0s=±l 

m iujQ 4- t{n, s) + z'A 

iujQ — e(n, s) + iA' m 



(30) 



with the free Landau level eigenvalue-spectrum 



Wo+e 2 (n,s) + m 2 = u% + \eH\(2n + 1 -s) + m 2 (31) 



(21) is A30h restricted to the LLL with n = 0, s = 1. 



Other matrix models with cross LL interactions are also 
possible, although suppressed by the large \eH\ gap be- 
tween LL. 



ERGODIC REGIME WITH STRONG 
MAGNETISM. 

For strong QED magnetic fields with a running cou- 
pling g(eH) that is not so small, the semi-classical con- 
figurations such as instant oris (vacuum) or calorons (ther- 
mal) are not suppressed. Although 4-dimensional, these 
configurations may considerably interact with the dimen- 
sionally reduced quark fields ^2 and localize them in near 
zero modes. A recent discussion of these near zero modes 
for a single instanton configuration can be found in |17j . 
Let N be the number of these near-zero modes in the 
transverse hyperplane to the magnetic field. Since the 
chiralities in 4 dimensions are commensurate with the 
chiralities in 2 dimensions, the pertinent chiral RMM for 
these near-zero modes is then standard 



*1* 4 ) = ^ I \eH\V*y/l- U */\eH\ (28) 



which is ( 23 ) for a symmetric box. The critical temper- 



ature is field dependent 



T,= 



VW\ 



(29) 



Ultra-strong magnetism implies tightly paired chiral 



quarks in the vacuum. The RMM results (23) and (29) 



Z LLL {m,T,A) = 

dAe -(A7A 2 )Tr(AtA) ^ 



m iujq + iA 

iojo + iA^ m 



(32) 

with A an Ax A matrix and A rj 1/p a typical QCD scale 
set by the instanton size p for instance. The gaussian 



weighting in ( 32 ) induces a 4- Fermi interaction 

2 



A- 



(*2* 2 y 



for the dimensionally reduced quarks. Note that the level 
spacing of the St 2 is 



AA 



A 
N 



(33) 



instead of (22 1. Again, by the Banks-Cashcr formula we 



expect the chiral condensate in a symmetric box to scale 
as 



7Tp(0) 



IT 2(j> 

ul AA 



I off I N 



A y/V 4 



^ (34) 



with n4 = N 2 /Vi « 1/fm 4 as typically the instanton- 
anti-instanton density in the QCD vacuum. 

The chiral condensate follows through similar reason- 
ing at the saddle point for large TV and at finite temper- 
ature 



^ T * 4 > = 



2N 



PL Z A 



^v^^ 



(35) 



N//3L Z is about the number of instantons and anti- 
instantons in f3L z regarded as a slice of V4 = (3L x L y L z . 
This is roughly \/n 4 since the instanton density 114 varies 
slowly with temperature except for T « T c where a re- 
arrangement into pairs is expected. Thus, the spacing 
( 33 1 between the eigenvalues is about AA ~ 1J^/V± in the 
presence of a strong magnetic field instead of AA sa I/V4 
in the vacuum. 

The chiral condensate is seen to grow linearly with the 
magnetic field. Unlike the precedent matrix model, the 
transition temperature is independent of the magnetic 
field 



T c = A/tt 



(36) 



Therefore, the slope of the chiral condensate for strong 
magnetism is 



*4*4 



VH4 



\eH\ y/1 (T/T c y 



(37) 



In the RMM, the chiral condensate at zero temperature 
is fixed by similar arguments. As a result, the relative 
change of the chiral condensate with magnetism is 



AE 



* 4 *4 



*4*4 



M y/1^ (T/T c f 
\/n,, 



(38) 



with n 4 sa 1/fm and T c about the critical temperature 
for chiral symmetry restoration without magnetism. 



We note that the effects of the higher LL on the quark 
condensate can be readily assessed in the generalized ran- 
dom matrix model ( p0| ) in the large <f> (ultra-strong mag- 
netism) and large N (strong magnetism) limit using the 
saddle point approximation. Indeed, in the latter limit 
and retaining the LLL n = 0, s = 1 along with the next 
to lowest LLL n = 0, s = —1 and n = 1, s = 1 yield 



N 



2\eH\ 

A 



1 - cj2/A 2 + 2Jl - 2\eH l/A 2 - lu 2 /A 2 



(39) 



The next to LLL contribution drops for strong fields or 
I off I > A 2 /2 as its corresponding saddle point moves off 
the real axis. Typically A w 0.5 GeV (see below) so that 
the decoupling takes place for \eH\ > 1/3 GeV 2 which is 
the onset of the ergodic regime. Thus our reduction to 
the LLL. 



COMPARISON WITH LATTICE SIMULATIONS 

For intermediate values of \eH\ < 1 GeV 2 our results 
for strong magnetism using random matrix theory and 
zero temperature are overall consistent with those re- 
ported by the recent lattice simulations in [11]. The de- 
viations at low values of \eH\ < 1/3 GeV 2 occur in the 
diffusive regime and are determined by the pion modes 
as we detailed earlier through the quark return probabil- 
ity. For ultra-strong magnetism our results appear to be 
consistent with the lattice results reported in [15] . 

At finite temperature, our results for strong magnetism 
using random matrix theory show a field independent 
critical temperature (36). The lattice results reported 



in |llj show a critical temperature weakly dependent on 
the magnetic field. The critical temperature appears to 
decrease by about 10% across the critical temperature. 
An anti-catalysis with a substantial decrease of the chiral 
condensate was reported to take place close to the critical 
temperature. 

In a recent investigation |12j the breaking of Z^ c sym- 
metry at high temperature with the appearance of a finite 
Polyakov holonomy was suggested as the mechanism for 
the appearance of the anti-catalysis phenomenon in the 
lattice data. The argument is that the finite temporal 
holonomy A4 = mp^ T with ^4 = (0, ±2/3) (mod 2) for 
N c = 3 or Lpi = 0, 1 (mod 2) for N c = 2 gets shifted 
by the magnetic field through sea effects. The result is 
a decrease of the chiral condensate around the critical 
temperature. In other words while the valence contri- 
bution is diamagnetic, the sea contribution through the 
Polyakov holonomy appears to be paramagnetic. 



Indeed, a small shift in the trivial Polyakov holonomy 
by the applied magnetic field can cause the chiral con- 
densate to decrease substantially near the critical tem- 
perature. The effect of the Polyakov holonomy is to alter 
the twisting of the temporal fermionic boundary condi- 
tion and therefore an up or down shift in the Matsubara 
frequencies. In the LLL the Polyakov holonomy is readily 
inserted through 



Z>LLL(m,T,A 4 



dAe -WA 2 )Tr(AtA) 



xdct 



i{u> 



rn 



iAt 



i(u 



m 



iA 



(40) 



instead of (32). The immediate effect of the Polyakov 



holonomy is to change the onset of the critical tempera- 
ture to T c — (A/7r)/|l + (fi\ (mod 2) in agreement with 
the result in [24] . 

The magnetic shift of the effective potential for the 
trivial Polyakov holonomy maybe estimated. At high 
temperature and for massless quarks, the shift in the 
pressure due to magnetism is [55] (see also section 6) 



An 



H 



l\eH\Hn' ^ 



\eH 



(41) 



with x = A r c/127r 2 . In the presence of a finite holonomy, 
(41 ) shifts through cj — > Wq+Aj (mod 2) with additional 



contributions that are not important for our arguments 
(see in general [23 ]). As a result, the Polyakov potential 
around the trivial holonomy (mod 2) reads 



AO, 



1 



? A 2 



— A 4 \eH\ 2 



(42) 



where we have added the leading contribution for zero 
magnetism. Here m 2 E is the electric screening mass. 
Without the electromagnetic shift, the minimum of (42) 



is around the trivial holonomy ip± = A^/uiq = (mod 2). 
The magnetic contribution causes this minimum to shift 



MH) « - X 



\eH\ 



(m E uo) 



(43) 



where the corrections are 0(e 3 ). We note that around the 
critical temperature in QCD with Nc = 3, the electric 
mass m E « ttT and x = l/^ 2 . So for T c ps 0.160GeV 
the Polyakov holonomy is shifted by about 10% for 
\eH\a (7rT c ) 2 «(0.5GeV) 2 . 



The shift (43) is rooted in the thermal cutoff in (41) 



through our perturbative estimate. Let us assume that 
generically the trivial holonomy shifts in the presence of 
magnetism with (fi(H) ps 0(e 2 ) as also suggested by cur- 
rent lattica data [T5]. While such a shift is small it can 



dramatically affect the chiral condensate near the crit- 



ical temperature. Indeed, inserting (43) in the random 



matrix model yields a chiral condensate 



* 4 *4 



N 

1 



2\eH\ 

A 



(44) 



]T Vl-(T/T c ) 2 (s + m (H)y 



=±1 



The sum over s = ± signs in ( 45 ) is important to un- 
derstand in our simplified random matrix theory. The 
chiral condensate in the saddle point approximation re- 
ceives contribution from all Matsubara frequencies u n = 
(2n + l)nT requiring in principle a matrix model with all 
frequencies [^ . Such a matrix model can be written 
and is explicitly periodic in 994 of period 2. However, near 
the critical points only the 2 lowest Matubara frequencies 
are dominant, i.e. luq = — lj_i (high temperature reduc- 
tion). In the absence of holonomies the random matrix 
model does not distinguish between these 2 frequencies 
near the critical points, thus the simplification to uiq. For 
a finite Polyakov holonomy, the 2 Matsubara frequencies 
are distinguishable, thus the sum over the signs. 

In the random matrix model the reduction of the chiral 
condensate takes place when 



Tg(H) 

T c (0) 



1 



1 - stp 4 (H) 



(45) 



where the subscript refers to the critical point. The w_ 



contribution in ( 45 ) drops first for a critical temperature 
T c (H)/T c (0) w 1 - \ifi(H)\. With the earlier estimates 
for \eH\ ps (0.5 GeV) 2 the change in the critical temper- 
ature T C (H) with the magnetic field is about 10% lower 
than without, in qualitative agreement with the lattice 



data [TT] . As a result ( 45 1 reduces at the critical point 
to 



*i*< 



N \eH\ 



-> 



(3L Z A 



Vftffl(l-W)l) (46) 



which is small and vanishes continuously with the in- 



crease of magnetism. We note that (45) implies a critical 



H C (T) through the perturbative estimate 



eH c (T) 



(m E uo) c 



1- 



T 



1/2 



(47) 



for the restoration of chiral symmetry at fixed magnetism 
and finite temperature for N c = 3. This is in qualitative 
agreement with the reported lattice data. The case of 
N c — 2 will be adressed in a sequel. 

As a final note, we would like to indicate that there 
is yet another way to test the nature of magnetism on 



the lattice and that is by twisting the fermionic bound- 
ary condition along L z with a Bohm-Aharonov flux A z . 
Specifically 



^> 2 {x Ql z + L z ) = -e l2 ^^ 2 {x Ql z) (48) 

This amounts to a bulk Abelian potential [201 EHJ &M 



(»7 • V E + A zlz ) * 4 = A[AJtf 4 



(49) 



with a constant A z — 2 r KLp z /L z . Recall that a constant 
A z is physical on a torus. Note that for (49) the quark 



virtualities A [A] are periodic in ip z of period 1. The 
twisted boundary condition ( 48 ) on the LLL in the ran- 



dom matrix model is achieved by using 



det 



m i(uj + A) + A z + iA 

iA+ + i(ui + At) — A z to 



(50) 



in (40 1. The chiral condensate follows through similar 



reasoning at the saddle point for large N with 



'* ' ' = W, 2 x^Vi-^/A 2 -K + A) 2 M 2 



tfj* 



(51) 
instead of (35). At zero temperature, the change in the 



chiral condensate is 



AS 



* *z 



\/u A 



(52) 



It would be interesting to carry numerical simulations 
of spectra and derive the chiral condensate from the per- 
tinent spectral fluctuations using directly the finite size 
random matrix construction with magnetism. Sea effects 
can also be analyzed this way through relevant rescaling 
limits [19 , although the random matrix model is totally 
dominated by the LLL with or without the Polyakov 
holonomy in the ergodic regime, as we have discussed. 
The effects of light quark masses and Np > 1 will be 
addressed elsewhere. 



Z LL (m,T,(t))= /"dBe-^ TrBtB det 



B + to iL 
iV Bt + m 

(53) 



with 



((wo + At) -ie(n,s))l T 



a diagonal matrix labeling the Landau levels shifted by 
the lowest Matsubara. The dominant contribution stems 
from the restriction to the LLL with l i and corresponds 
to the reduction *&± to ^ ■ The addition of the spatial 
twist is straightforward. The corresponding effective ac- 
tion for the LLL in 2-dimcnsions is 



9. 



LL 



N 



TrB f B- Alndet 



B + to iL 

iV Bt + to 



(54) 



(54) restricted to the LLL describes the dimensionally 



reduced effective action that describes the quark spec- 
trum near zero virtuality for strong fields in the ergodic 
regime. For instance 



*J*: 



lim lim 

m— fO L z — foe 




B + to ill 

iV B f + to 



(55) 



The averaging is carried using the measure in ( 53 1 . The 
content of this dimensionally reduced effective action 
along with the role of a spatial twist and a real quark 
chemical potential will be discussed elsewhere. 

In what follows we will develop a more phenomenolog- 
ically inspired random matrix effective action that inter- 
polates between the diffusive and ergodic regime through 
a reduction of the NJL model to constant modes [3] (and 
references therein). The results are inspired chiral RMM 
effective actions with magnetism. The purpose of this 
construction is to show the potential relationship between 
the random matrix theory discussed above as suggested 
by the mesoscopic analysis and the body of work on mag- 
netism using NJL inspired models. 

To analyze the chiral Euclidean effective action or pres- 
sure we recall that the Gibbs energy for finite T and fi 
of chiral RMM is that of a two-level system for a parti- 
cle and antiparticle for a single Matsubara frequency and 
Np = l [S[TS] 



RMM EFFECTIVE ACTION WITH 

MAGNETISM. 

The random matrix effective action with magnetism 



follows from ( 30 1 through standard fermionization of the 



determinant and then bosonization by the use of the 
Hubbard-Stratonovitch transformation. The result for 



the partition function ( 30 ) is 



n 4 



2 



^ ln(( f T- M ) 2 + (7rT) 2 )(( f 7 + Ai ) 2 + ( 7 rT) 2 )) 

(56) 

with a playing the role of the constituent quark mass and 
E 2 the strength of the chiral Gaussian noise. It ties to 
the chiral condensate in the vacuum through 



(*+* 4 ) = sV, 



(57) 



at the extremum. Again, the density n 4 = N/V4 refers in 
general to the number N of zero modes in the 4-volumc 
which is commensurate with the density of instantons 
plus antiinstantons in the vacuum. 

AB fluxes act as Abelian U(l) vector potentials and 



are readily implemented in ( 56 ) through the substitution 



(i^ 2 ^(iH + A±) 2 + A 2 z (58) 

We recall that on a torus constant gauge fields are gauge 



invariant. For /1 = A4 = T = (56) simplifies 



IW(A a ) = i£V-5i]n(<7 



2 ' A 2 ) 2 



(59) 



with a chiral condensate at the minimum 



*i* 4 



£ 2 <7* = SV2n 4 - (A Z S) 2 (60) 



that decreases with the twisting of the boundary condi- 
tion as we noted earlier. 

The virtuality levels of constituent quarks in a fixed 
external magnetic field along the z-direction are fixed 
by (14 1. The vacuum energy follows by summing over 



the negative energy states through standard procedure. 
Specifically, the negative of the energy per unit V4 or 
pressure reads 



n H (o,o) 



2 



\eH 



(2tt 



OO 

»E 

n=0 



dk z (2-5 n0 )\- tl (0,k z ) 

(61) 

where we have added the external magnetic field energy 
density contribution. The sum is UV divergent and re- 
quires renormalization. Most importantly, we note that 
in the vacuum state all the states are paired in spin with 
the exception of the state n = and a = 1. The disor- 
dered smearing of this near-threshold state is responsi- 
ble for the spontaneous breaking of chiral symmetry by 
strong magnetism as we showed above. 

The leading order in 1/N C effects to the free con- 
stituent quark loop is the dominant contribution whereby 
the constituent quark mass is generated through weaker 
residual and attractive interactions say of the instan- 
ton type. For a fixed external magnetic field, the con- 
stituent quark loop contributes to the magnetic permit- 
tivity n(H). For a fixed cut-off Ajjv and a strong mag- 
netic field or A uv > y/\eH\ > (a + m) [22] 



fiff(0,0) = 



H 2 



MH) 



ff 2 

2 



%e' 



hi 



A 2 
eH 



(62) 




FIG. 1: Magnetic contributions to the constitutive quark 
loop. See text. 



for Np — 1. For several quark flavors e 2 — > ^W F e g with 
Np = 3, e u = 2e/3,ed — — e/3,e s = e/3. Since the loop 
contribution stems from the UV sector, the constituent 
quarks are for all purpose massless, making the logarith- 
mic contribution unique. The plus sign of the vacuum 
contribution reflects on the diamagnetism and hence the 
stability of the vacuum state under the switching of a 
constant H. All the remaining contributions from the 
loop are analytic in 0((eH) 2 ) and go into the charge and 
H renormalization. The combination eH is renormaliza- 
tion group invariant. We set the renormalization so that 
for Afjy = \eH\ (62) is canonical [2"2"] . 



The result ( 62 ) is readily understood from the fermion 
loop diagram of Fig. [l] with the fermionic lines corre- 
sponding to a constitutive fermion moving in the exter- 
nal magnetic field. The free loop with 2 H-insertions is of 
the form H 2 /p with p a typical loop momentum. Clearly 
the integral diverges both in the ultraviolet and infrared. 
The cutoff in the ultraviolet is Ajjv while in the infrared 
it is a + m for free constitutive fermions and 



(a + m) 2 -> \eH\ + (a + mf 



(63) 



for constitutive fermions in a magnetic field. For Auv > 
\/\eH\ > (cr+m) (62) follows, while for A uv > (a+m) > 
\Z\eH\ we have 



fiff(0,0) 



H 2 



We- 



ill 



A 2 



(a- 



(64) 



For \J\eH\ > (cr+m) we note that the higher H-insertions 
in Fig. [T] are increasingly infrared sensitive making them 
all of order |eiJ| 2 . They renormalize the classical field 
contribution 122 . In dense matter near the critical 



points, the UV cutoff in (64) is traded through 



A 2 



\ojq + i(i\ 



(65) 



with ujq = 7rT, leading to 



9 




to B (T,n) 



(66) 
in agreement with detailed limits derived in [22j. (66) 
permeates all constitutive quark calculations in the con- 
text of the magnetic catalysis of chiral symmetry break- 
ing in QCD both in vacuum and matter [9] (and refer- 
ences therein). 

Our starting point for a RMM analysis of the phase 
diagram at finite H in matter with twisted boundary 
conditions is fi = CIm + CIh- With only H present, the 



minimum dfl/dcr* = using ( 56 ) and ( 66 1 yields 



cr* =E v / 2n 4 + X^ 2 



(67) 



which is to be compared with ( 60 ) . This yields a change 
in the chiral condensate 



AS = 




- 1 



(68) 



which is quadratic for low \eH\ and linear for interme- 
diately large \eH\. We recall that for ultra-strong mag- 
netism at zero temperature the ratio asymptotes \eH\ 3 / 2 
as both suggested by RMM and perturbative QCD. 



The asymptotic of ( 68 1 



AS 




is smaller than ( 38 1 by the factor 



y/Nc/24/ir « 1/10 



(69) 



(70) 



The random matrix model is not restricted by the 1-loop 
result as it involves the derealization of the dimension- 
ally reduced chiral quarks through quasi-2-dimensional 
zero modes and is more in line with the currently reported 
lattice results. The results of the random matrix analysis 
are reproduced by setting x — ► 2 in fl H . This effective 
action analysis without boundary twists whether tempo- 
ral (Polyakov) or spatial (BA), emphasizes the underly- 
ing diamagnetic character of the random matrix analysis 
above when restricted to the LLL. 



CONCLUSIONS AND PROSPECTS 

We have presented a chiral RMM for the analysis of 
the recent lattice QCD simulations in the presence of a 
strong QED magnetic field. The results for intermediate 



values of \eH\ > l/p 2 ~ 0.3 GeV 2 are readily understood 
through the infrared branch of the LLL which causes the 
quarks to dimensionally reduce from 4 to 2 dimensions. 
Quarks trapped in 4 dimensional LLL are more prone to 
break spontaneously chiral symmetry through random 
disorder by gluonic configurations whether perturbative 
or non-perturbative. This is the essence of the magnetic 
catalysis discussed in [9j (and references therein). Note 
that our reduction of the spectrum is in quark virtualities 
not quark energies and therefore from 4+1 to 2 + 1 rather 
than 3 + 1 to 1 + 1. 

In 2 dimensions the interactions due to the semiclas- 
sical vacuum configurations, i.e. instantons or calorons 
cause them to accumulate at zero virtuality. The ac- 
cumulation is characterized by a level spacing of order 
1/y/Vi which is intermediate between I/V4 in the vac- 
uum and 1/ y/Vi in a free box. The change in the chiral 
condensate is shown to increase linearly with \eH\ (catal- 
ysis) with a slope l/y'nj « 1 fm with 114 the instanton 
vacuum density. These results are consistent with the 
current lattice data at zero temperature. For very weak 
\eH\ the gapped charged pions dominate the contribu- 
tion in the diffusive regime in line with the leading result 
from chiral perturbation theory. At ultra-strong mag- 
netic fields \eH\ > 10/ p 2 w 3 GeV 2 the RMM suggests 
an increase in the chiral condensate as \eH\ 3 ^ 2 in agree- 
ment with perturbative QCD estimates for ultra-strong 
fields using a Bethe-Salpeter analysis [13"] . 

At finite but high temperature the emergence of a 
Polyakov holonomy does not affect the decoupling of the 
LLL. We have shown that a small shift in the trivial 
part of this holonomy allows for a rapid decrease (anti- 
catalysis) of the chiral condensate near the critical tem- 
perature as suggested recently by lattice simulations [H] . 
This small change in the Polyakov holonomy is quadratic 
in the magnetic field. It reflects on the back-reaction 
of the quarks on the effective potential of the temporal 
holonomies near the trivial solution. 

We have explicitly shown how the random matrix effec- 
tive action emerges from our analysis in the dimension- 
ally reduced limit. Also, we have provided a specific con- 
struction of a number of random matrix inspired models 
to allow for a simple comparison with current constituent 
quark models with magnetism such as the NJL models. 
A more detailed analysis of the phase diagrams emerging 
from these effective actions will be detailed elsewhere. 

In random matrix theory the spectral distributions in 
the microscopic limit capture more information on the 
subtleties of chiral symmetry breaking such as the role of 
quark representations, quark masses and magnetism and 
twists. 

Of particular interest is Nc = 2 with an external mag- 
netic field which upsets time-reversal symmetry, a big 
deal for this representation. Indeed, for Nc — 2 the chiral 
disorder involves not only the pions (diffusons) but also 
the baryons (cooperons) both of which are degenerate fol- 
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lowing the extra Pauli-Gursey symmetry. The cooperon 
is sensitive to the breaking of time-reversal symmetry be- 
cause of its charge [2D] ■ 

Finally, it would be interesting to revisit the role of 
the number of flavors, quark masses and the spatial and 
temporal twists on the results we have derived for a more 
critical comparison with current lattice data. Also the ef- 
fect of the chemical potential on the LLL analysis should 
shed more light on the role of magnetism on the QCD 
phase diagram. Some of these issues will be addressed 
next. 
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